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Abstract. We study the number of solutions of the general semigroup equation 
in one variable, X a — X 13 , as well as of the system of equations X 2 — X 1 Y 2 = 
Y, XY = YX in H lT n , the wreath product of an arbitrary finite group H with the 
full transformation semigroup T n on n letters. For these solution numbers, we provide 
explicit exact formulae, as well as asymptotic estimates. Our results concerning the 
first mentioned problem generalize earlier results by Harris and Schoenfeld (J. Combin. 
| Theory Ser. A 3 (1967), 122-135) on the number of idempotents in T n , and a partial 

result of Dress and the second author (Adv. in Math. 129 (1997), 188-221). Among 
f-H , the asymptotic tools employed are Hayman's method for the estimation of coefficients 

■ of analytic functions and the Poisson summation formula. 
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^ ■ 1. Introduction 

00 

(Sj ! For a finite group G and for a positive integer d, denote by sc(d) the number of 

subgroups in G of index d. It is well known that the number | Hom(G, S n ) \ of G-actions 
^ ■ on an n-set satisfies the identity 

m ' 

§: f^|Hom(G,S n )|^ = exp Ij^^-zA , \G\ = m; (1.1) 

n=0 "■ \d\m j 

cf., for instance, ^21 Prop. 1] or 6, Prop. 1]. Formula (jl.lj) . which exhibits the ex- 
ponential generating function of the sequence | Hom(G, ^n)! <^> <^ 

rather simple type of 

entire function, can be made the starting point for the asymptotic enumeration of finite 
^ [ group actions; cf. [13, Sec. 1] or [Tlj . 

Of course, after choosing a finite presentation for G, \ Hom(G, S n ) \ can be interpreted 
as the number of solutions in S n of a certain system of equations (corresponding to the 
set of relations in this presentation). For instance, denoting by ^#(m) the number of 
solutions of the equation X m = 1 in a finite group H, the case G = C m of (jl.lj) yields 
the formula 

oo „ / \ 

5>„(m)^ = ex P XXA* , (1-2) 

n=0 ' \d\m J 



first proved in ^j. 
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Moreover, let G and H be finite groups. Then it is shown (among other things) in 
[IS] that 

£ | Hom(G, H I S n )\^ = exp ( £ HtMffi^ ) , (1.3) 

n=0 ' \d|m " / 

where 

So(d):= E |Hom(C/,#)|, 

(G:U)=d 

this formula leading to corresponding asymptotic estimates for | Hom(G, H\S n ) \ by the 
methods developed in [Tl] . A noteworthy special case, originally due to Chigira [3], is 
the formula 

n=0 \d|m / 

In this case, the main result of ^1] provides an asymptotic expansion for Lms n { m )) 
whose main term yields the asymptotic formula 

L mSn ( m )~K(\H\nf-^ex V (-(l-±) n + ^^^n d A , (1.5) 

V d\m ' ' J 

where 

y /m, m odd, 

K - \ . — / |/./ \nt\-n \ 

-y/m exp — — — — , m even. 



The present paper arose out of an attempt to establish analogues of ()l.lj) - (jl.5j) for 
the number of solutions of equations (or systems of equations) in the transformation 
semigroups T n , and, more generally, in wreath products H\T ni where if is a fixed finite 
monoid. 

Our first result, Theorem ^ computes the exponential generating function for the 
number s(n) of solutions in H \T n of the most general semigroup equation in one 
variable, that is, 

X a = X p , 0<a<p, 
in the case where if is a finite group. We find that 

n>0 ' ^' 1 7 |/3-a ' 

where the sequence of functions {A a (z)} a >o is recursively defined via 

A (z) =z, 

(1.7) 

A a (z) = 2;exp(A Q _i(2)), a > 1. 

Moreover, Theorem Q also computes in closed terms certain refinements of the series 
J2 n >o s ( n ) zn / n - Thus, it extends a previous result by Dress and the second author 
Prop. 2] from transformation semigroups to wreath products of groups with transforma- 
tion semigroups. The special case of (jl.6|) where H = 1 answers a question of Stanley's 
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posed in fU\ Ex. 6.9]. We list several other interesting specializations as examples after 
the proof of Theorem 

Combining with Hayman's machinery ^T] and some detailed asymptotic cal- 
culations then leads to our first main result, namely an asymptotic formula for s(n) 
in terms of the real root of a certain transcendental equation; cf. Theorem |3] Formu- 
lae (jl.6|) and (j3.12|) should be seen as semigroup analogues of Formulae (jl.4j) and (jl.5j) 
in the group case. We conclude the first part of the paper by applying this general 
formula in the special case a = 1; see the example after the proof of Theorem In 
particular, we recover the result by Harris and Schoenfeld |9, Eq. (29)] concerning the 
asymptotics of the number of solutions of X 2 = X in the transformation semigroup T n . 

For systems of equations, the situation appears to be considerably more involved. 
The (larger) second part of our paper is concerned exclusively with an analysis of the 
particularly simple system of equations 

X 2 = X, (1.8) 

Y 2 = Y, (1.9) 

XY = YX. (1.10) 

We begin by combinatorially characterising the solutions of (jl.8|) - (|l.l(J|) in H I T n , 
where H is a finite group, and then use this characterisation to explicitly compute the 
exponential generating function for the number a(n) of these solutions. The result is 



^aH^ = exp(^|ifr 1 cH^), (1.11) 

where c(n) is equal to 



n\ \ * — ' n\ 

n>0 \ra>0 



C (n) = — r n-r-s-t s n-r~s+l. (1-12) 

^— ' r! s\ t\ (n — r — s — t)\ 

0<r,s,t<n v ' 

cf. Corollary [7| 

The second main result of our paper is the determination of the asymptotics for a(n). 
In Theorem ITUl we find that 

a{n) = \H\ n ~ l c{n) (l + O (-)) , as n -> oo, (1.13) 



where 



n 



n 1/2 3/2 
U \ Tri Sri 



c(n)~e s «^ 2 ' — as n — > oo, (1.141 

\s n J n 

with (r n ,s n , t n ) the unique solution for large n of the system of transcendental equa- 
tions (|7.5|) - (|7.7|) . The functions r n , s n , and t n are only implicitly defined, and their 
asymptotic description naturally involves iterated logarithms as gauge functions, and, 
hence, are difficult to approximate. We prove ()1.13|) by determining the asymptotics 
of the triple sum ()1.12|) (leading to ()1.14j0 . and then combining the latter result and 
the generating function identity (jl.llj) with Bender's transfer method pQ. The hardest 
part of this argument concerns the asymptotic estimation of c(n) , the reason being that 
classical analytic methods based on generating functions appear not to be applicable. 
Instead, we work directly from the triple sum (|1.12j) . determine the location of those 
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(r, s, t) which provide the main contribution to the sum, and compute its value by 
applying the Poisson summation formula; see the proof of Theorem |H1 

Comparison of our result for the asymptotics of a(n), the number of solutions of 
(I1.8|) - (jl.l0|) in H I T n , with our result for the asymptotics of the number of solutions 
of X 2 = X in H I T n reveals that the effect of commutativity is, asymptotically, almost 
negligible, as it can be seen only in the third term of the asymptotic expansion; compare 
Remark El in Section El 

Our paper is organised as follows. In Sectional we consider the exact enumeration 
of solutions of the equation X a = X 13 in HlT n , where if is a finite group. We solve the 
problem in Theorem Q by providing a compact formula for the exponential generating 
function of (a refined version of) these solution numbers. The subsequent section, 
SectionEl is devoted to determining the asymptotics of these solution numbers, which we 
accomplish in Theorem El The remaining sections deal with the problem of enumerating 
solutions (X,Y) of the system of equations ()1.8J) - (jl.lOJ) . In Section HJ we concentrate 
on solutions in T n , and solve the problem of exact enumeration of these solutions in 
Theorems 0] and Corollary El by studying carefully the combinatorial structure of such 
solutions. Then, in Section El we extend these findings to the enumeration of solutions 
of this system of equations in H lT n , where H is a finite group; see Theorems El and 
Corollary [7| Sections EHU deal with the asymptotics of these solution numbers. The 
"opening" section, Section El provides an overview of the results, the main result being 
Theorems El and El and outlines their proofs. Missing details of these proofs are then 
supplied in Sections [7| and El 

2. The equation X a = X 13 in H l T n 

For a non-negative integer n let T n denote the full transformation semigroup (or sym- 
metric semigroup) on n letters, that is, the set of all maps on the standard n-set 
[n] = {1,2, ... ,n} with multiplication given by (ri ■ r 2 )(j) := r 2 (ri(j)). For a (finite) 
semigroup H with unit element 1 we form the semigroup 

HlT n = {(f,r):feH^, r G T n }, 

with multiplication given by 

(fi,n) ■ (/ 2 ,t 2 ) = (/,ri • r 2 ), 

f(j) = /i(j')/ 2 (ria)), je[n). 

Under this multiplication rule, H I T n is a semigroup with unit element (l,id), called 
the wreath product of H with T n , the map 1 being identically 1. In what follows, H 
will always be taken as a finite group. 

Given r G T n , the set [n] decomposes into non-empty r-invariant subsets Sj, [n] = 
Si U 5*2 U ■ ■ ■ 5/;, such that each Sj is indecomposable, that is, Sj itself is not the union 
of two disjoint non-empty r-invariant subsets. This decomposition of [n] is uniquely 
determined by these requirements, and is nothing but the partition of [n] given by the 
vertex sets of the connected components of r viewed as a directed graph on [n]. The 
equivalence relation inducing this decomposition of [n] is given by 

ji ~ 32 there exist integers k, I > such that r k (ji) = T l (j 2 ). 

We call Si, S 2 , ■ ■ ■ , Sk the components of r, the numbers \Si\, \S 2 \, ■ ■ ■ , \Sk\ are the de- 
composition numbers, and k is the decomposition length of r. Given a pair of sets A, M, 
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where A C N and M C No, a map r : [n] — > [n] is termed (A, M)- admissible if the 
decomposition numbers of r are in A and the decomposition length is in M. 

Given a finite group H, integers a and j3 with < a < /3, and sets A C N and 
M C No, we write Sjf(n) for the number of (A, M)-admissible solutions of 

X a = X? (2.1) 

in H I T n . 

Our first result computes the exponential generating function for the numbers s^(n) 
in the case where if is a finite group. It extends a previous result by Dress and 
the second author [HJ Prop. 2] from T n to wreath products H I T n . For a power series 
f( z ) = J2k>o f kzk an d a se t K of non- negative integers we define the restriction {f{z))x 
to be the series Ylk&K fk z>t - 

Theorem 1. Let H be a finite group. Then we have 

= e "(M ^ 7-M^)((A a (l#k)) 7 ) A ), (2-2) 

n>0 ' W \ 7 | /3 _ a / 

where Ln{ m ) is the number of solutions in H of the equation x m = 1, ejvf(z) := 
(exp(z))jvf; ° nc ^ ^« as (|1.7|) . 

Proof. A pair (/, r) E H \ T n satisfies (|2.ip if and only if 

(i) r satisfies (12. lj) in T n , and 

(ii) for all j G [n] the equation /KO')) = Iltl f{r v {j)) holds in H. 
Since if is a group, (ii) is equivalent to 

(ii') nt^(^j)) = 1 for all J G [4 
A solution t of ()2.1|) in T n is the disjoint union of oriented cycles of lengths dividing 
(3 — a, with ingoing trees of height at most a (possibly 0) attached to each cycle point; 
and conversely, each map r G T n of this form is a solution of (j2.1j) . See Sec. 3.2] for 
more details. 

We have to determine, for each such r, how many / G satisfy (ii'). Fix r. If 
j G [n] occurs in one of the trees of r, then f(j) can be chosen arbitrarily in H, since 
such a vertex does not occur in (ii'). If c is a cycle of length 7 in r, and if xi, x 2 , ■ ■ ■ , x 7 
denote the images of the vertices of c under /, then Equation (ii') for these vertices is 
equivalent to the condition that 

(x 1 x 2 ---x^~ a)h = 1. (2.3) 

The problem of solving (J2.3)) can be decomposed into first choosing x G H such that 
x {P-a)h = i ) anc l then solving the equation x\x 2 ■ ■ ■x 1 = x. Hence, Equation (|2.3J) has 
precisely 

solutions. It follows that each solution r of (|2.ip in T n can be paired with 

ni^'-MM) (2-4) 

if 
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maps / G such that (/, r) satisfies (|2.1|) in HlT n , where K runs over the components 
of r, and c(K) is the cycle of K. Hence, writing |r| := n if r G T n , we have 

£•?<») ir^S II w^-'-dSS). < 2 - 5 > 

n>0 T K component of r 

where the sum is over all r whose number of components is in M, where the sizes of 
the components are in A, and where the lengths of the cycle in the components divides 
(3 — a. 

The best way to get a manageable expression for the exponential generating func- 
tion of the numbers sjjf(n) in ()2.5|) is to think in species-theoretic terms (cf. [2] for 
an in-depth exposition of "species theory"). The transformations r which appear in 
the sum on the right-hand side of ()2.5|) belong to the (species-theoretic) composition 
Sets M (Cycles /3 _ Q! (Trees< a )), where Sets A ,/ denotes the species of sets with cardinal- 
ity in M, Cycles /3 _ a denotes the species of cycles of length dividing (3 — a, and where 
Trees< Q denotes the species of trees of height < a. The species Setsj\/(Cycles /3 __ a 
(Trees< Q )) contains in fact a larger class of objects, since the condition that the sizes of 
components should be in A is not yet taken into account. So, let (Cycles /3 _ Q ,(Trees< Q ,)) 
denote the subspecies of Cycles /3 _ Q (Trees< Q ,) which consists only of those objects hav- 
ing a size which is contained in A. Then to compute the right-hand side of (J2.5)) 
amounts to computing the (weighted) exponential generating function for the species 
Sets M ((Cycle 

s /3-a( Trees <a))A), where the weight of an object from that species is 

given by (|2.4j) . 

Clearly, the exponential generating function for Setsj^ is ejw(x), the one for 
Cycles /3 _ ce is J2y\/3- a zl lit an d the one for Trees< a is A a (z) with A Q (z) as defined 
in (jl.7|) ; the last assertion following inductively from the exponential principle (cf. e.g. 
[21 Theorem 1.4.2(a) with F = E, the species of sets, so that F(x) = exp(x)]). Hence, 
by the substitution theorem for species (cf. e.g. (21 Theorem 1.4.2(a)]), the exponential 
generating function for Cycles /3 _ Q ,(Trees< Q ,) is 

However, we need a weighted enumeration of objects K from that species with the 
weight 

depending on the size \K\ of the object and the size |c(if)| of the (unique) cycle of K. 
Hence, we have to modify ()2.6|) appropriately. Indeed, the weighted generating function 
for Cycles /3 _ Q (Trees< ct ) with weight ()2.7j) is 

1 ^ /0 _ a , (A a (\H\z)r 

W\ 1^ lh \—) — - — • 

Now, this series has to be restricted to powers z m with m G A to obtain the (weighted) 
generating function for (Cycles /3 _ Q ,(Trees< Q )) A , which, upon using the substitution 
theorem once more, must then be substituted into e^iz), the generating function for 
SetSjvf. This will give the weighted generating function (|2.5j) for the species Sets M 
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((Cycles /3 _ Q (Trees< a ))A) that we are looking for. The result is exactly the right-hand 
side of □ 

Examples. For A = N, M = No, and H — 1, Formula (|2.2j) specializes to the formula 
[HI Sect. 3.3.15, Ex. 3.3.31] of Goulden and Jackson, which in turn answers a question 
posed in [21]. However, the right-hand side of (|2.2j) is readily calculated for every choice 
of A, M, and H; for instance, we find for a = 1 that 

E*A(»)^ = e*( £ 7 -W-M^)^(A-W7l#k)Y (2-8) 

n>0 ' 71/3-1 ' 

As a further illustration, let us take a closer look at the case a — 1, /3 = 2 of idempotents 
in if 2 T n . Setting sf (n) = f/ n (A, M, H), we find from (J22J) that 

^ 17„(A, M,H)— = e M (z e A ^(\H\z)) . (2.9) 



n>0 



Taking M = 2N and A = 2N - 1, the right-hand side of flHH|) becomes 

cosh(zcosh(|/f |,)) = i + l(|/f| 2 ) 2 + ^(i^k) 4 + ^(im 6 + ^(\m 8 +■■■, 

while choosing M = 2N + 1 and A = 2N yields 

sinh(z sinh(|.ff|2:)) 

2 ,. , .„ 4 ,, , 4 126,, __. , fi 3368, . , s 95770, . , in 

= 2? (I^k) + l*0 4 + -gi-(l^k) 6 + -^-(l#k) 8 + "i^-d^i^ + ■ ■ • • 

As a final example, let us take M = 2No and 

A = A = { n g N : n = 1 (3) and n > 4}. 

By ()2.9j) . the exponential generating function of the number U n (A , 2No, H) of idempo- 
tent elements in H I T n with an even number of components all of whose lengths are 
congruent to 1 mod 3 and > 4 is given by 

2>(A.,2N.,Jf)£r = cosh E^i^r 

n=0 \re=l v 

„ 560, lrrl , 8 9240, , n 124124,, „. , 14 672672000 ,, „. , 18 

= i + ^r(l^k) 8 + -^(Itf I*) 11 + -^j-(l^k) 14 + — —i — (l# k) 16 + • • • . 

3. ASYMPTOTICS OF THE SOLUTION NUMBERS FOR X a = X 13 IN H I T n 

Here we shall combine Theorem with analytic machinery developed by Hayman to 
analyse the asymptotic behaviour of the function s(n) := s^°(n) counting the number 
of the solutions of (J2HJ) in H I T n . Setting A = N and M = N in (|2~2|l . the exponential 
generating function of s(n) is seen to be equal to 

*(*) := ^ s(n) J = exp (±- £ ^ Lh{ ^) (A a (\H\ z) y) . 
The version of Hayman's result best suited for our purposes is the following. 
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Theorem 2. (Hayman) If f(z) = ^2 n>0 fnZ n is admissible in \z\ < R in the sense of 
[TT] . then 

f n ~ , n as n — > oo, (3.1) 

where r n is the unique solution for large n of the equation a(r) = n in (Rq,R), with 
a{r) = rf'(r)/ f(r), b(r) = ra'(r), and a suitably chosen constant Rq > 0. 

This is PH Cor. II, p. 69] (see also [THJ Cor. 12.6]). 

It follows immediately from Hayman's definition of admissibility that zf(z) is admis- 
sible in \z\ < R provided f(z) is admissible in this range. Combining this observation 
with the facts that e z is admissible with R = oo (see [TH1 Example 12.7]) and that 
admissibility is preserved under exponentiation ( Theorem VI]; ^Sl Theorem 12.8]), 
we see that the functions A a (z) defined in (|1.7|) are admissible, for each a > 1, in the 
whole complex plane . The Corollary to fTJ Theorem IX] (see also |18| Theorem 12.8]) 
combined with the above exponentiation rule now shows that ^{z) itself is admissible 
in the complex plane. 

From now on we shall assume that a > 1. Setting f(z) = ^(z) in the context of 
Theorem El we find that 

a(r) = r £ tff (5z2) (A^Hlr))^ 1 A' a (\H\r) } (3.2) 

j\/3-a 

and that 
6(r) = a(r) + \H\r 2 £ t H (^) 

x h-l)(A Q (|i/|r)) 7 ~ 2 (A^(|i7|r)) 2 +(A Q (|i7|r))^ 1 A' c :(| J ff|r)]. (3.3) 

As in Theorem El let r n be the real root of the equation a(r) = n. Next, we shall 
obtain an asymptotic estimate for r n . Clearly, r n —>■ oo as n —>■ oo; hence, for r = r n , 
the dominating term on the right-hand side of (J3.2j) is the term with 7 = /3 — a, and 
thus 

n^r n {A a {\H\r n )Y~ a ' l A' a {\H\r n ) asn^oo. (3.4) 
Taking logarithms, this gives 

logn ~ logr n + (J3-a- 1) log A a (\H\r n ) + log A' a (\H\r n ). 
Moreover, it follows by induction on a that 

A' a (\H\r n ) = A a (\H\r n )g(n), 
where exp(g(n)) = 0(A a (\H\r n )). This leads to 

logn ~ logr n + ((3 -a) log A a (\H\r n ). 
Since a > 0, we may also neglect the term logr n , thus obtaining 

logA a (\H\r n ) ~ — — log(n). 

p — a 

Substituting A a (\H\r n ) = \H\r n exp (A a _i(|if |r n )) , we get 

log(|#|r n ) + A^fliflTvO ~ ^5 log(n). 

p — a 
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Neglecting the first term on the left-hand side and taking logarithms gives 

logA a _i(|#|r n ) ~ log (— --— log(n)). 

For convenience, let us write \og^ m \z) for the "m-th iterate" of the logarithm, that is, 
\og^ m \z) := \og(\og^ m ~ l \z)) and \og^\z) := z. Then, iterating the above procedure, 
we find that 

Ao(|#|r) ~ log^ 1 ) (-A_log(n)), 

and thus, finally, 

r " ~ w\ log(a_1) log(n) ) ' a - L (3 ' 5) 

Note that the factor l/(/3 — a) in (|3.5|) can be dropped for a > 2. 

Putting r = r n in (13. 3|) . and using ()3.4|) as well as the fact that the dominating term 
on the right-hand side of (j3.3J) is the summand for 7 = (3 — a, we find that, as n — > 00, 

A'«(i#K) , A£(|tf|r n ) 



6(r n ) ~ n 



l + |#|r„ [((3-a-l) 



By definition of A a (2;) we have 

A'Jz) 



A a (z) 
This gives 



A Q (\H\r n ) A' a (\H\r r , 
z^ + A'^z), a>l. 



(3.6) 



A a -i(z) as z — > 00. (3.7) 



A a (z) A Q _!(2; 
Iterating (|3.7j) and using the fact that 

A'^z) ~ Ax(z) as 2; — > 00, 

we find that 

^^-A^A^-.-A^) as z^oo. (3.8) 

Again, by the definition of A a , 

K(z) = -z~ 2 A a (z) + z-^'Jz) + A^A'^z) + A'^Kiz). 
Dividing both sides of this equation by A' a (z), and using (|3.7jl . we get 

Iterating the last equation, using (|3.8|) and ignoring neglectable terms, we obtain 
A"(z) 

Ai(^)A a (z)---A«_i(z) as 2^00. (3.9) 



A' (2) 



Setting 2; = |-H"|r„ in ()3.8|) and ()3.9|h and using the resulting relation in ()3.6|) . as well 
as (|3.5jl . we obtain 

b(r n ) ~ (^-ajnlogl"- 1 ' (_L_l og (n)) A^tflrgA^lEK) • • • A^O^K). (3.10) 
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Using (|3.2j) and the fact that a(r n ) = n allows us to rewrite ^(r n ) as 

*(r n ) = exp 



nA a (\H\ 




(J3-a)\H\ 


r n K(\H\ 





Inserting this equation and (|3.1U|) in Hayman's formula 1)3.1 Jl finally gives the following 
result. 

Theorem 3. Let H be a finite group, and let a and (5 be integers satisfying < a < (3. 
Then, as n — > oo, the number s(n) of solutions of the equation X a = X 13 in H \T n is 
asymptotically equal to 

s(n) 1 = 

CyXS - a)A 1 {\H\r n )A 2 (\H\r n ) ■ ■ ■ A a ^(\H\r n ) \og^ log(n)) 

(( A a ([g|r n ) \ 

(3-12) 

7|/3— a 

where r n is the unique solution for large n of the equation 

r E ^(^)(A Q (|i/|r)) 7 - 1 A' a (l^)=n. 

y\/3-a 

Remark 1. For a = 0, the function s(n) in Theorem El counts the number of solutions 
of the equation X 13 = 1 in the group HlS n . The asymptotics for this case, and, indeed, 
for the more general case of the function Hom(G, H I S n ) with finite groups G and H, 
can be found in [T^l Theorem 4]; cf. also [TBI Eq. (30)] or [T71 Eq. (49)]. 

Example. For a = 1, the statement of Theorem |3] simplifies, and we find that the 
number s(n) of solutions in H I T n of the equation X 13 = X is asymptotically equal to 



where r n is the unique solution for large n of 

r ^ ^(£_I)(|#| r )7-l( 1+ |fr| r ) e 7|H|r = n _ 
7 |/3-l 

If we further specialize to the case where (3 — 2, the last formula simplifies further to 



n 



s(n) ~ 



— ■ 



er n J Vlogn 
where r ra is the unique solution for large n of 

r(l + r)e r = n \H\. 
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In the case that H = 1, this result has been previously obtained by Harris and Schoen- 
feld; cf. |3 Eq. (29)]. For later use we observe that asymptotically, as n — > oo, we 
have 

(log log n \ 
— , 
logn J 

and, thus, 

/ \ i / \ , n v l0 ^log(logn) n /n(log(logn)) 2 \ 

sin) = n log in)—n log (logn)— n+2 h YO 5 . 3.13) 

logn logn ^ (logn) 2 J 

4. The system of equations X 2 = X, Y 2 = Y, XY = YX in T n 

In this section we consider the system of equations (|1.8|) - ()1.1U|) . and we ask for the 
number of solutions (X, Y) of this system in T 2 . 

Setting a = 1 and (3 = 2 in the proof of Theorem [H we see that a solution X of 
(jl.8J) in T n , when regarded as a directed graph on [n], is the disjoint union of loops 
with arrows pointing to the vertex of the loop; see Figure ^ for such a component. Of 
course, the same applies to any solution Y of (jl.9j) . 




Figure 1. 



Let (X, Y) be a solution of the system (jl.8j) - (jl.l0j) . We regard such a pair (X, Y) as 
a directed graph on [n], where the edges can occur in two colours: the colour red, say, 
associated with the edges of the graph corresponding to X, the other colour, blue say, 
corresponding to the edges of the graph corresponding to Y. 

As a first step, we are going to enumerate the solutions (X, Y) in T 2 whose directed 
graph is connected, by which we mean that any two vertices can be reached from each 
other by travelling along edges of either colour and regardless of the direction of the 
edges. A typical example is displayed in Figure where red edges are symbolized by 
fully drawn arrows, and blue edges are symbolized by dashed arrows. Let us write c(n) 
for this number. 

Theorem 4. For n > 1, the number c(n) of connected solutions (X,Y) e {T n ) 2 of the 
system (I1.8j) - (jl.l0|) is given by 

C(n) = Y , , ,, Ul n-r-s-t s n-r-s+l ^ 

K J ^ r«t n-r-8- 

0<r,s,t<n v > 

where 0° has to be interpreted as 1, and where ml is interpreted as oo for negative values 
of m, so that the sum in (14.1)1 is actually restricted tor + s + t<n. 

Proof. As we already observed, we may regard any solution (X, Y) of the system (jl.8j) - 
(jl.lUJl as a graph with red and blue edges, where the 'X-components" (connected 
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components consisting of red arrows only) are loops with arrows pointing to the vertex 
of the loop, which are "connected" by blue arrows. 

Let us consider such an 'X-component," and denote its X-fixed point by v; see 
Figure ^ We ask ourselves: which are the possible images of v under Y (that is, to 
which vertex could v be connected by a blue edge)? We consider first the possibility 
that v is mapped under Y to a vertex, u say, different from v and in the same X- 
component, that is, to a vertex u with X(u) = v, see Figure 121 In this case we would 
have Y(X(v)) = u and X(Y(v)) = v , a contradiction to (jl.lOj) . 




Figure 2. 



A second possibility is that v is mapped under Y to a vertex w in another X- 
component. This vertex w may be a fixed point under X or not. If w is not a fixed 
point, then X(w) = u for some u 7^ w. Then we encounter the situation illustrated 
in Figure El In that case, we have Y(X(v)) = w and X(Y(v))) — u - - again a 
contradiction to (|1.1U|) . 




Figure 3. 



Hence, the only possibility is that fixed points under X are mapped to each other 
under Y . 

Now we consider one of the vertices in an X-component, which is not fixed under 
X, the vertex w, say. Can this vertex be mapped under Y to another vertex that is 
not fixed under X, to the vertex t, say (see Figure H|? Let the X-fixed point of the 
component containing w be denoted by v, and let the X-fixed point of the component 
containing t be denoted by u. Since t may lie in the same X-component as w, we may 
also have u = v (in that case we would have to replace the picture in Figure |U by the 
appropriately modified picture). In either case, because of X(Y(w)) = Y(X(w)), it 
follows that u = Y(v). In summary, whenever vertices which are not X-fixed points 
are mapped to each other under Y, then the X-fixed points in the corresponding X- 
components must be mapped to each other. 

Consequently, if (X, Y) is connected, then there is exactly one X-fixed point, t> say, 
which is also a fixed point of Y, and all other fixed points of X are mapped to vq under 
Y. Thus, the typical form of a connected solution (X, Y) is as shown in Figure El There 
are several X-components, K , Ki, . . . , K r , say, whose X-fixed points (in the figure 
these are v , v±, . . . , v r ) are mapped under Y to the same X-fixed point v . In the 
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Figure 



X-component K , the component containing vq, there may be some vertices different 
from vq which are mapped to themselves under Y. (For example, in Figure^ the vertex 
U\ is such a vertex.) Let the number of these vertices be s — 1 with some s > 1, and let 
the number of the remaining vertices in K Q (including t> ) be t + 1. Out of the latter, 
the t vertices different from vq are mapped under Y to one of the s fixed points under 
Y (including v ). This is indicated in Figure by the dashed lines going out of w 2 and 
m 3 . Finally, the vertices in the other X-components K\, . . . , K r which are not fixed 
points under X can be freely mapped under Y to the s fixed points under Y in Kq. 
Examples in Figure El are the vertices u^, uq, or u§. (The reader must imagine that also 
the unfinished edges in Figure El lead to one of the s fixed points under Y in Kq.) 




Figure 5. 



Now let r, s, and t be fixed. How many such connected solutions of (jl.8)) - (jl.ltJ|) are 
there (that is, configurations with n vertices such as in Figure EJ), which have exactly 
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r+1 X-components, which have s+t vertices in the X-component containing the vertex 
which is at the same time a fixed point for X and Y, and which have s fixed points of 
Y in that X-component (including the vertex which is fixed under X and Y)? 

Firstly, out of n possible vertices, we have to choose the r fixed points under X which 
are not at the same time fixed under Y (in Figure El these are the vertices V\, . . . , v r ), 
we have to choose the s fixed points under Y in the X-component K (in Figure 03 these 
are Vq and, e.g., ui), we have to choose the t other vertices in K Q (in Figure El these 
are, e.g., «2 and U3), and we have to choose the remaining (n — r — s — t) vertices. The 
number of these choices equals the multinomial coefficient 



r! s\ t\ (n — r — s — t)\ 

Next, having made this choice, we must mark vq, the vertex which is to be the fixed 
point under both X and Y, out of the s fixed points under Y in K . Clearly, this allows 
for s possible choices. Having done all that, the only freedom left for X- and Y-images 
of vertices is to determine the X- and Y-images of the (n — r — s — t) "remaining" 
vertices and the Y-images of the t vertices in K different from the Y-fixed points. The 
candidates for the X-images of the (n — r — s — t) "remaining" vertices are v\, . . . , v r . 
Hence, there are exactly r n ~ r ~ s ~ t choices. The candidates for the Y-images of the 
(n — r — s — t) "remaining" vertices and the t vertices in Kq different from the Y-fixed 
points are the s fixed points under Y in Kq. Hence, there are exactly s n ~ r ~ s choices. 

In total, we obtain 

. g . j.n—r—a—t . g n-r-s 

r! s\t\ (n — r — s — 1)\ 

choices. If we sum this term over all possible values of r, s, and t, then we obtain the 
expression in the statement of the theorem. □ 

Corollary 5. Letb(n) denote the number of all solutions (X, Y) e (T n ) 2 , of the system 
(OJ) - (IPTll) . Then 

2>(n)— = exp X)c(n)— , 

n>0 \n>0 ' / 

where c(n) is the number determined in Theorem® 

Proof. This follows immediately from the exponential principle (cf. e.g. [3 Theorem 
1.4.2(a) with F = E, the species of sets, so that F(x) = exp(x)]) and Theorem □ 

5. The system of equations X 2 = X, Y 2 = Y, XY = YX in H I T n 

Here, we shall generalize the result of the previous section from the symmetric semigroup 
T n to the wreath product HlT n , where H is a finite group; cf. Section[2]for the definition 
of the wreath product. 

Our next theorem determines the number d(n) of connected solutions (X, Y) of the 
system of equations (jl.8j) - (jl.l0j) in H I T n . Being connected here means that, if X = 
(/, a) and Y = (g, r), then the pair (a, r) is connected in the sense of Section HJ As in 
Section |U we denote by c(n) the number of solutions in T n of the system (jl.8|) - (jl.l(J|) . 
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Theorem 6. The number of connected solutions (X, Y) G (H I T n ) 2 of the system 
(jl.8|) - (|l.l(Jj) is given by d(n) = \H\ n ~ 1 c{n) , where c(n) is the number determined in 
Theorem® 

Proof. Let (X, Y) G (HlT n ) 2 be a solution of the system ljn ty -(|TT0 |) . By the definition 
of the wreath product, this is equivalent to 

a 2 = a, t 2 = r, or = ra, (5.1) 

and 

/(<) " /(*(<)) = /(»*), (5-2) 

^(i)-^(r(i)) = ^(i), (5.3) 

/(0 •»(*(<)) = »(<)• /MO), (5-4) 

the latter equations being required for all i G [n]. Equation (J5.2j) is equivalent to 
f(j) = 1 for j G Imcr. Similarly, Equation (J5.3|) is equivalent to g(j) = 1 for j G Imr. 
If i G Imcr, i = a(j) say, then (|5.4j) becomes 

m-g(vHj)))=g(i)-f(T(o-(j))). 

Since a 2 = a, ar = ra, and in view of our previous observations, this equation becomes 
1 ■ g(i) = g(i) ■ 1, a tautology. Similarly, if i G Imr, Equation ()5.4|) is automatically 
satisfied by our previous findings. 

We assign arbitrary values to f(i) for i G Imr\Imo" and to g(i) for i G Imcx\Imr. 
If we write r = | Imr\ Im a\ and s — 1 = | Im a\ Imr| (which is in accordance with the 
notation in the proof of Theorem then these are r + (s — 1) assignments, for which 
we have \H\ r+s ~ 1 possibilities. We shall see that no contradiction arises if this is done. 

Now let i £ Ima U Imr. We claim that f(i) can be chosen freely, and determines 
g(i) uniquely. Indeed, in Equation (|5.4|) the values g(a{i)) and f{j{i)) have been 
already assigned, and therefore f(i) determines g(i). Since every element in Imcr is 
in fact a c-fixed point, and similarly for r, and since we are considering a connected 
solution, the intersection Imcr fl Imr consists of exactly one element (compare the 
proof of Theorem HJ). Hence, we have | Imcr U Imr| = r + s, and, consequently, there 
are |if | n ~ r ~ s possibilities for assigning values of / and g for % ^ Imcr U Imr. 

In total, given a connected pair (cr, r) satisfying (j5.1J) . there are 

|^-|r+s— 1 |^|n— r— s | | n — 1 

possibilities to assign values for / and g. The claim of the theorem now follows imme- 
diately. □ 

Corollary 7. Let a(n) denote the number of all solutions (X, Y ) G {H I T n ) 2 , of the 
system flUEj) - ((HDD ■ Then 

£a(n)- = ex P n>>rMn)- , 

n>0 \n>0 '/ 

where c(n) is the number determined in Theorem® 

Proof. This follows immediately from the exponential principle (cf. e.g. [21 Theorem 
1.4.2(a) with F = E, the species of sets, so that F(x) = exp(x)]) and Theorem |3 □ 
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6. ASYMPTOTICS OF THE SOLUTION NUMBERS FOR THE SYSTEM X 2 = X, Y 2 = Y, 

XY = YX IN H I T n 

In the remainder of this paper our aim is to establish an asymptotic estimate for the 
number a(n) of solutions (X,Y) G (H I T n ) 2 of the system of equations (jl.8|) - (jl.l0|) . 
Note first that the exponential generating function of the sequence c(n) (and hence also 
that of a(n)) diverges outside the origin. Indeed, by considering the summand in ()4.1|) 
with r = s = |_?V(21ogn)J and t = 0, and applying Stirling's formula, one sees that 

. c(n) 1 , 

log — > -72 log 

n\ 2 

Furthermore, attempts to use other types of generating functions, for instance 

n>0 

appear to lead to intractable analytic functions. Thus, complex analytic tools like 
Hayman's in Theorem El do not seem to be applicable in this context. Instead, we will 
work directly from the sum formula for c(n) given in Theorem EJ and, in combination 
with a theorem of Bender, we shall obtain a compact asymptotic expression for a(n); 
see Theorem ITUl below. 

The crucial point is the proof of the following result, the main steps of which will be 
described below, leaving technical details for later sections. 

Theorem 8. We have 

, x n 1/2 3/2 

C(n) ~ e ^{n-r n + \)~Zn+2r n +2s n+ t n +l I ^ \ r " g " ^ 

\s n J n 

„ e ^(n-r n +i)-3n+2r n +2 Sn+ t n +l ( ™\ " as 71 — > OO, (6.2) 

\s n J 41og 2 n' V ; 

where (r n , s n , t n ) is the unique solution for large n of the system of equations (J7.5|) - (j7.7|) . 

Outline of proof . Given a large positive integer n, consider the function 

Y(n + 1) 7" n ~ r ~ s ~* s n—r—s+l 

F(r ' S ' t] = T(r + 1) r(s + l)T(t+l)T(n-r-s-t + l) (6 ' 3) 

for r, s, t e Q n , where 

tt = Ur,s,t) E R 3 : r > 0, s > 0, t > 0, r + s + t < n\. 

In Lemmas fT4l and ITBI we demonstrate that F(r, s, t), which is nothing but an extension 
to real numbers of the summand in (j4.1j) . attains its global maximum at a unique point 
(r n ,s n ,t n ) in the interior of Q n . This maximum point (r n ,s n ,t n ) is the solution of 
the system of transcendental equations (|7.5|) - ()7.7|) . and Lemma HU also provides the 
estimates 

n 

r n ~ s n ~ — : , as well as t n ~ 21ogn. (6.4) 

2 logn 

Next, in Lemma |2*HI we show that contributions to the sum (j4.1j) . which in our present 
notation is 

0<r,s,t<n 
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which come from outside of a certain box Q n containing (r n , s n , t n ) in its centre, are of 
the order of magnitude 

o(F(r n ,s n ,t n )n/ log 2 n). (6.5) 
On the other hand, the summands of the "smaller" sum of summands coming from Q n , 

E F(r,s,t), 

(r,s,t)eQ n 

can be approximated in terms of more tractable functions, so that a consequence of the 
Poisson summation formula, given in Lemma l2"Ul can be applied to obtain the compact 
estimate 

p( (2ir) 3/2 r n s ny /t n ~ 
n 

for the latter sum; cf. Theorem [IM] By (|6.4j) . this term dominates the one in (J6.5j) . The 
final step consists then in finding an asymptotic expression for F(r n ,s n ,t n ), which we 
do in Proposition [213 If everything is combined, the first line of formula (j6.2j) results. 
Obviously, the second line is a simple consequence, if one uses ()6.4|) . □ 

Remark 2. As was said at the end of the last proof, the somewhat more explicit expres- 
sion in the second line in ()6.2|) follows immediately from the expression in the first line, 
given the asymptotic information on r n , s n , and t n provided by ()6.4|) . We have stated 
the first expression, since the second one has considerably worse approximation be- 
haviour, due to the fact that the quotient of r n by its asymptotic equivalent n/(21ogn) 
converges to 1 only at a sub- logarithmic rate, with a similar statement applying to s n . 
More details can be found in Remark HI and Equations (|7.9|) - (|7.11|) . 

Remark 3. We expect the relative error of the asymptotic formulae in Theorem [H] to 
be of order 0(l/n). A proof of this would require considerable extra work: at various 
places in the proofs of the Lemmas in Sections [7| and |S1 the effects of higher terms in 
the occurring expansions would have to be taken into account, leading to a massive 
calculational effort. 

The translation of (|6.2|) into an asymptotic estimate for the function a(n) is made 
possible by the main results of Bender [JJ (see also [TBI Theorem 7.3 with F(x, y) = 
exp(y)], as well as [23 121) • 

Lemma 9. // the power series A(x) = J2 n>l a n x n and B(x) = ^ n>1 b n x n satisfy 

(i) B(x) = exp(A(x)), 

(ii) a n -i = o(a n ) as n — ► oo, 

n-l 

Q'k&n— k 



E 



k=l 

then 



n-l 



0(1) 



n-l j 



We are now in the position to state the second main result of this paper, and to 
describe the main steps of its proof, again leaving technical details for later sections. 
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Theorem 10. Let H be a finite group, and let a(n) denote the number of solutions in 
H \ T n of the system of equations (jl.8j) - (jl.K)j) . Then 

a(n) = \H\ n ~ 1 c(n) (l + O (-)) , asn^oo, 



with the asymptotics of c(n) being provided by Theorem® 

Outline of proof . We use Lemma El with a n = c(n)\H\ n ~ 1 jn\ and b n = a(n)/n\. In view 
of Corollary the corresponding series A(x) and B(x) satisfy (i). It remains to verify 
the conditions (ii) and (iii). 

We first consider (ii). Using the approximation ()6.2|) in Theorem |H1 we obtain 

a„_i n c(n — 1) 
a n \H\ c(n) 

1 / n — 1 / 1 , 

■ exp (n-1) - r n _i + - - 3(n - 1) + 2r n _i + 2s n _i + £ n _i + 1 



\H\ V s„_i \ 2 / 

+ (n — 1) log(n — 1) — (n — 1) log s n _i + log(n — 1) — 2 log log(n — 1) 

n ( 1\ 
— I n - r n + - ) + 3n - 2r n - 2s n - t n 



— nlogn + nlogs n + 2 log log n I . (6.6) 

We now appeal to the asymptotic information on r n ,s n ,t n provided by Lemmas fill 
and[T7J which we substitute in the expression (|6.6|) . The result is then simplified; for 
instance, one would estimate the term (n— l) 2 /s n _i —n 2 /s n , occurring in the exponent 
in (16. 6J) . as follows: 

(n-1) 2 n 2 i n - l f s n~ n2 Un~^ + o(^ 



2 1 

n fT^ 1 - 2n , . 
= — H h o (log n) 

■Sn— l^n Sn— 1 

= (2 logn + o(logra)) — (4 logn + o(logn)) +o(logn) 

= — 2 logra + o(logn). (6.7) 

After a tedious but straightforward computation, we obtain that 

~ — , as n — > oo, (6.8) 

a n n 

with an explicit constant C±, whose exact value is of no relevance here. This establishes 
(ii). 

In order to estimate the sum 

n-1 

occurring in (iii), we observe first that it suffices to consider the sum from k = 2 to 
k = [n/2\, since the sum is symmetric with respect to the substitution k — > (n — k), 
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and since the term for k = 1 is just 1. Now we divide the latter range into two parts, 
the interval Ii := {k : 2 < k < log log n} and the interval I2 '■= {log log n < k < n/2}. 
By Theorem |H1 and Lemma IT41 we have 

a m = exp (m log m + o{m log m)), as m -> 00. (6.10) 

Therefore, if k 6 Ji, then a fc < exp ((log log n) (log log log n)) <C yjn. This implies 



the last inequality being due to ()6.8|) . Therefore, the sum 



E 



is bounded. 

If k 6 then, using Theorem we have 

a k a n __ k = Iffp ^^'f ~ l#P ^p ($„(*;) - 2n + 2 - 2 log4 - ± log(27r)) , 
k\{n — k)\ 

with $ n (/c) the function defined in Lemma fTUl Hence, uniformly in k and n, it is true 
that 

a k a n _ k < C 2 |#| n exp ($„(&) - 2n) , 

for some constant C2. Now, by Lemma ITH| we have $«(&) < $n(2) for /c e /2- Using 
Theorem El again, and the definition of $ n (2), it follows that 

a k a n ^ k (n- 1)! exp($ n (2) - 2n) 

a n -i c{n - 1) 

< C 3 exp I ^ — - (^n - r n _ 2 - ^ + 2r n _ 2 + 2s n _ 2 + t n „ 2 

- (n - 2) log s n _ 2 + - log(n - 2) - 2 log log(n - 2) 



■Sn-l V 2 



n - ?V_i - - - 2r n _i - 2s n _i - 



+ (n - 1) log s n _i - i log(n - 1) + 2 log log(n - 1) 

for some constant C3. By a similar calculation as the one which showed that the 
expression (jfi.7|) is of the order of magnitude 1/n, we deduce that a k a n - k / a n -\ < C^/n, 
for some constant C4. Hence, the sum 



E 



, , a n-l 
keh 

is also bounded, and, thus, the complete sum (|6.9|) . as required. □ 

Remark 4. Combining Theorems IH1 and ITTfl we obtain an asymptotic formula for a(n) in 
terms of the transcendental functions r n , s n , t n . In view of the speculative Remark 01 
this formula should give an approximation with a relative error of 0(l/n). On the 
other hand, by inserting the asymptotic expansions (|7.9|) - (|7.11|) for r n , s n , t n into 
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Equation (JO}, one obtains an asymptotic expansion for the exponent of a(n) in terms 
of elementary functions. The first terms of the latter expansion are 

tXl log log TL 
2n log n — 2n log log n — 2 (log 2 + 1) n + 3 — : 
logn 

+ (3 log 2 + l)-^ + of " (1 ° glo y )4 ))^ (6.11) 
logra y (lognf J J 

The gauge functions occurring in the exponent of 1)6.11)) decay at a sub-logarithmic rate; 
hence, while providing a rough idea of the order of magnitude of a(n), Equation (J6.ll)) 
is not useful for numerical purposes (whereas (J6.1)) in combination with Theorem 
is). 

Remark 5. It is interesting to compare the asymptotics for the number of solutions of 
()1. 8)1 - 1)1.10)1 in H I T n on the one hand, and the square of the number of solutions of 
X 2 = X in the same sequence of semigroups, since it allows us to estimate to what extent 
the commutative law ()1.10)) restricts solutions of ()1.8)) and (J1.9)) . Using the expansions 
()6. 1 1(1 and (J3.13)) we find that the quotient of these two numbers is asymptotically 

a(n) ( nloglogn n /n(loglogn^ 4 
- exp — (21og2)n h (3 log 2 — 1)- h O 1 



s(n) 2 \ \ogn logra \ (logn) 2 

hence the commutative law effects only a restriction of exponential order on a sample 
space growing roughly like n 2n . 

7. Auxiliary lemmas I: existence, uniqueness, and asymptotic properties 

of (r n , s n , t n ) 

The main results of this section are Lemma ^1 where we prove uniqueness and 
determine the location of the global maximum (r n ,s n ,t n ) for the function F(r,s,t) 
defined in (16. 3)1 . Lemma ITTJ1 in which we prove that the auxiliary function <& n (k) defined 
there has a unique minimum and no maxima, and the asymptotic informations on 
r n ,s n ,t n gathered in ()7.8)1 and Lemma fTTl All other results in this section are of a 
preliminary nature, leading to the proofs of the aforementioned results. 

Our first two lemmas deal with location and asymptotics of the maxima for an aux- 
iliary function f(r,s). 

Lemma 11. Given n G N, consider the function 

T(n + l) r n-r-s s n-r-s+l 



f{r, s) 



r(r + 1) r(s + 1) T(n - r - s + 1) 

for (r, s) G Vt n , where Q n := {(r, s) GM 2 :r>0, s>0, r + s < n}. Then, for 
sufficiently large n, f(r,s) attains its maximum in the interior ofQ n . 

Proof. Since Q n is compact and f(r, s) is continuous on Q n , f attains a global maximum 
on Q n . We will show that 

n n 



f{r,s) < f 



2 log n 2 log n 
for (r, s) G dQ n and n large, which implies our claim. 
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Application of Stirling's formula shows that 

/ ( — r— — , — ; ] = exp (2nlogn + (9(n log logn) ], as n — > oo. (7.1) 

\ 2 log n 2 log n J \ / 

On the other hand, we have 

T(n + l)s n ~ s+1 _ n+1 n 



/(0, s) < ^ ' ' -r < 2V +i « / 



r(s + l)r(n-s+ 1) V 21 ^^ 21ogn 

and an analogous estimate holds for f(r, 0). If r + s = n, then 

r , s T(n + 1) s , / n n 



T(r + l)r(s+l) ' \2 logn 2 logn 



□ 



Lemma 12. Let (£ n ,??n) € fi n 6e sucn £/iet£ f(£ n ,Vn) is maximal on Q n . Then 

Kn ~ ??n ~ 771 ' OS 72 ► OO. 

2 logn 

Proof. Since / attains a maximum at (£ n , 7/ n ), and as (£ n , r] n ) dVL n by Lemma fTTI the 
point (£ n ,Vn) satisfies the equation (v/)(6n> ^n) = 0, which is equivalent to the system 
of equations 

- log£ n + - — y — — - \ogr) n - ?/>(£„ + 1) + ip(n -in- TJn) = 0, (7.2) 
-log£„ + - — — — ^ + 1 _ logn» - ip(r] n + 1) + ip(n - £ w - rjn) = °) ( 7 - 3 ) 

Vn 

where ip(x) denotes the digamma function. 

Suppose first that £ n ^> nj logn. Then the term (n — £ n — r] n )/i n m f)7.2|) is <C logn, 
whereas the remaining terms give a negative contribution of absolute value 3> logn, 
which is seen upon recalling the asymptotic behaviour 



^(x) = \ogx-^-j^ + o(J-^, asx^oo, (7.4) 

of the digamma function (cf. 7, 1.18(7)]). But this contradicts (|7.2|) . Hence, £ n < 
n/ logn, and a similar argument based on Equation ()7.H|) shows that r\ n < nj logn. 

Now suppose that £ n <C n/ logn. Then the term (n — £ n — r) n )/£ n in (|7.2j) is ^> logn, 
whereas the remaining terms are at most O(logn), again contradicting (J7.2|) . Hence we 
have £ n x n/logn, and the analogous argument based on Equation ()7.3|) shows that 
r] n x nj logn. 

Now let £ n = a n (nj logn) with bounded ot n . Then another glance at Equation ()7.2|) 
gives 

log n 

= 2 logn + O(loglogn), as n — > oo, 

that is, a n = | + o(l), and thus £ n ~ n/(21ogn). 

The same argument with Equation (J7.3|) yields n„ ~ n/ (2 logn). □ 

Our next two results are analogues of Lemmas fTTI and IT21 for the function F(r, s,t). 
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Lemma 13. Given n e N, consider the function 

T(n + l) r n-r-s-t s n-r-s+l 



F(r,s,t) 



T(r + 1) r(s + 1) T(t + 1) T(n - r - s - t + 1) 

for r, s, i G O nj where 

ft = {(r, s, t) E M 3 : r > 0, s > 0, t > 0, r + s + 1 < n}. 

Then, for sufficiently large n, F(r, s, t) attains its global maximum at points in the 
interior offt n . 

Proof. We argue in a similar fashion as in the proof of Lemma ^2 in view of (j7.1|) we 

have 

— , — ,0 =/ 771 ' TTi = exp 2n log n + 0(n log log n) . 

2 log n 2 log n J \2 log n 2 log n J \ / 

On the other hand, we have 

V ' ' ' - T(s + l)T(t + l)T(n - s - t + I) - \2\ogn 21ogn' 

and an analogous estimate holds for F(r, 0,t). If r + s + t = n, then 

p ^ s ^> r(r + l)r(s + l)r(* + l) - V21ogn'21ogn'7 ' 



Finally, we have 



F(r, s, 0) = F(r, s, 1). 

n — r — s 



In Lemma ITT1 we proved that the left-hand side attains its maximum in the interior of 
Q n . Furthermore, by Lemma IT21 any such point {iniVn) satisfies £ n ~ r\ n ~ n/(21ogn). 
Hence, for large n, 

F(r,s,0) < F(UVn,0) = _ f n _ F(UVnA) < F(^ Vn ,l). 

Thus, we have shown that, for n large, F(r,s,t) cannot attain its maximum on the 
boundary of Q n . This establishes the claim. □ 

Lemma 14. Let (r n , s n , t n ) G fi n be such that F(r n , s n , t n ) is maximal on Q n . Then 

- log r n - log s n + - — — — — — — - ip(r n + 1) + 4>(n - r n - s n - t n + 1) = 0, (7.5) 

- log r n - log s n + - — — — — - ip(s n + 1) + ip{n - r n - s n - t n + 1) = 0, (7.6) 

S n 

-\ogr n - ip{t n + 1) + i/j(n - r n - s n - t n + 1) = 0. (7.7) 
Moreover, we have 

n 

r n ~ s n ~ — and t n ~21ogn. (7.8) 

2 logn 
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Proof. Since F attains a maximum at (r n ,s n ,t n ), and since (r n ,s n ,t n ) ^ dQ n by 
Lemma Uni the point (r n ,s n ,t n ) satisfies the equation (\/F)(r n , s n ,t n ) = 0, which is 
equivalent to the system (|7.5JI - (|7.7|) . 

To obtain the claimed asymptotic estimates, we first assume that t n x n. Then 
Equation (|7.7|) implies that r n = 0(1). On the other hand, the term {n—r n — s n —t n )/r n 
in (J7.5p must be < logn, because by f|7.4|) all the other terms in ()7.5j) are. Thus, the 
left-hand side of (17. 7J) is in fact < —logn + O(loglogn), a contradiction. Hence, we 
have t n <C n. 

Now the arguments from the second half of the proof of Lemma IT21 can be used again, 
where instead of Equations (|7.2j) and (|7.3j) one considers Equations (|7.5|) and (|7.6|) . to 
show that r n ~ s n ~ n/(21ogn). 

Assuming that t n = 0(1), the left-hand side in (j7.7J) would be asymptotically equal 
to log logn + o(loglogn), a contradiction. Hence t n ^$> 1. If we now substitute all this 
information in f)7.7j) . then we obtain logt n = log(2 logn) +o(l), whence t n ~ 2 logn. □ 

Remark 6. Further terms in the asymptotic expansions of r n , s n , t n can be obtained by 
an iterative procedure. Having already obtained the first few terms in the expansions, 
one adds, for each of r n , s n , t n , a further indeterminate. This "extended" expansion 
is then substituted into Equations (|7.5j) - (J7.7|) . Subsequently, asymptotic expansions 
are determined for the left-hand sides of these equations. Inspection of leading terms 
then yields a system of linear equations for the three indeterminates, which is solved. 
Thereafter, the next loop of the procedure can be started. Clearly, it is not advisable 
to do this by hand. However, with the help of Salvy's Maple programme gdev jl 9| |2"U] , 
these calculations are conveniently performed. The first terms, obtained in this way, 
turn out to be as follows: 

1 n 3nlog(logn) (3 log 2 — 2) n 9 n (log(logn)) 2 
rn ~21og^ 4 (logn) 2 + 4: (logn) 2 + 8 (logn) 3 

18 log 2 - 21) nlog(logn) (9(log2) 2 - 21 log 2 + 8) n 



+ - ^ , + 



8 (logn) 3 8 (logn^ 3 

27 n(log(logn)) 3 (162 log 2 - 243) n (log(logn)) 5 



+ T77 " "„ + 



16 (logn) 4 32 (logn) 4 

(81 (log 2) 2 - 243 log 2 + 135) log(logn) 



+ 



16 (logn) 

(54(log2) 3 - 243(log2) 2 + 2701og2-80) n ( n (log(logn))^ } 

+ ^ 7, vi +U \ 7] ^5 h U- y J 

(logn) \ (logn) 



1 n 3nlog(logn) 31og2 — 2 n 9 n (log(logn)) 2 

21og^ 4 (logn) 2 + ^ (log(n)) 2 + 8 (logn) 3 

18 log 2 - 21) nlog(logn) (9(log2) 2 - 21 log 2 + 8) n 



+ ^ , + 



8 (logn) 3 8 (logn^ 3 



■fi y/ i <- 

(logn) 
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, 01 01 n n 01 o , , 91og(bgn) , (91og2-4) 

t n = 2 logn — 31og(logn) — 31og2 — 1 + -- 



2 logn 2 log n 

+ 27 (log(logn)) 2 + (27 log 2 -39) log(logn) + Q ( 1 
8 (logn) 2 4 (logn) 2 \(logn)' 

Lemma 15. For n large enough, the function F(r,s,t) attains its global maximum on 
Q n at a unique point in the interior offl n . 



Proof. Let (r, s,t) be a point in the interior of Q n with r ~ s ~ n/(21ogn) and t ~ 
2 logn, and let e = (a,b,c) be a unit vector in IR 3 . We consider the function F along 
the line segment (r, s, t) + Ae, with A varying over a suitable interval containing 0. 
Differentiating twice with respect to A we find 

d 2 

— ( log F(r + Xa,s + \b, t + Ac)) 
= -na + b + c)^ a - v ^(n-r-s-t-(a + b + c)X) 

-2(a + b)^ b - I ^ w (n-r-s-(a + b)X + l) 

- a 2 ^ (1) (r + Aa + 1) - b 2 ^ 1] (s + Xb + 1) - c 2 ^ (1) (t + Ac + f) 

- (a + 6 + c) 2 ^ (1) (n-r-s-t- (a + 6 + c)A + l), (7.12) 

where ip^\x) is the first polygamma function. 

For A <C n/ log n the second and fourth term on the right-hand side of ()7.12j) have 
order of magnitude (logn) 2 /n, while all other terms are 0(logn/n), which is seen upon 
making again use of the asymptotic expansion for the first polygamma function (cf. 
1.16(9) and 1.18(9)]), 

^W(ac) = - + ( — ) , asx^oo. (7.13) 



X \x 2 



Hence, for large n and A -C n/ logn, we have 

d 2 



( log F(r + Aa, s + Xb, t + Ac)) < 0. 



dX 2 

This shows that along each line segment 

n 

[r,s,t) + Ae : A < 

logn 

the function F attains at most one maximum for n large enough. 

Now suppose that F attains its maximum at two different points (r 1; si,ti) and 
(r2,S2,t 2 ). Then, by Lemma ITU r\ ~ r 2 ~ S\ ~ s 2 ~ n/(2logn) and ti ~ t 2 ~ 2 logn. 
Considering the line segment through (r 1; Si,ti) and (r 2 ,S2,t 2 ) we obtain r x = r 2 , 
Si = s 2 , and ti = t 2 , contradicting our assumption that (r 1; Si,ti) 7^ (r 2 , s 2 ,t 2 ). □ 

In the sequel we write r' n for the derivative of r n with respect to n, regarding n as a 
real variable. (It can be checked that the proofs of Lemmas ITTHTol concerning existence 
and uniqueness of the solution (r n , s n ,t n ) for large n of the system of equations (|7.5j) - 
(17. 7|) remain valid for real n. Hence, r n , s n , and t n are well-defined for large real n.) 



EQUATIONS IN FINITE SEMIGROUPS 



25 



Similarly, the symbol r" n denotes the second derivative of r n with respect to n. We use 
the notation s' n , s", t' n , and t" n with analogous meaning. 

The next two lemmas, Lemma Hoi and IT%| provide asymptotic estimates for r' n ,s' n ,t' n 
and r^,s^,t", respectively These are needed in the proofs of Corollary IT71 and Lem- 
ma the latter two being crucial for the proof of Theorem ITU1 

Lemma 16. Let (r n ,s n ,t n ) e Q n be such that F(r n , s n ,t n ) is maximal on Q n . Then 

1 . , 2 



n °n 



Moreover, we have 



nr r 



2 logn 



and t' 



n 



n 



ns, 



2 log 2 n ' 



as n 



as n 



oo. 



oo. 



(7.14) 
(7.15) 



Proof. We consider the system of equations (|7.5j) - (|7.7|) for real n. Then, by the implicit 
function theorem, the system will have a (local) solution (r n ,s n ,t n ) as long as the 
determinant of the Jacobi matrix of the system is non-zero. Now, with E X ,E 2 ,E 3 
denoting respectively the left-hand side of (|7.5jl . (|7.fij) . and f!7.7|) . and writing r, s, i for 
r n , s n , respectively t n , the Jacobi matrix J of the system is 

fdE 1 dE l 8E X \ 



J 



dr 


ds 


dt 


dE 2 


dE 2 


dE 2 


dr 


ds 


dt 


dE 3 


dE 3 


dE 3 



Here. 



dE x 

dr 
dEi 

ds 
dE x 

~df 

dE 2 

dr 
dE 2 

ds 
dE 2 

dt 
dE 3 

dr 
dE 3 

ds 
dE 3 

dt 



\ dr ds 



"' + — - ^W(r + 1) - ^ {1) (n -r - s-t+1), 



1 

r 



ip^ 1 ' {n — r 



-- - ^ (1) (n - r- s- t + 1) 



r 
1 

r 



i, 

s 



r 

71 



n 
71 



1 1 

-V> (1) (^-^-s-t + l) 



+ 1) -V (1) (n-r- + 



tb^fn — r 

r 



s-t+1), 
s-t + 1). 



—ip^'in — r 
-ij {1) (t + l) - (n - r - s 



t + 1). 



Using the asymptotic information on r n , s n ,t n provided by Lemma ITU and the asymp- 
totic expansion (|7.13j) for the first polygamma function, it is straightforward to check 
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that det J ~ 2/ logn as n — > oo. This implies that det J must be non-zero for almost 
all n. Hence, again by the implicit function theorem, for all n > N (where iV is an 
appropriate real number) the solution (r n , s n ,t n ) of the system (|7.5j) - (|7.7j) exists and is 
differentiable, and the derivatives r' — r' n , s' — s' n , t' — t' n satisfy 

/- 



J 



- V W (n 

- ipW(n 



r - s - t + 1)\ 

r-s-t + 1) . (7.16) 

\ —ipw{n - r - s - t + 1) y 

Solving this system for r', s' and t', and using again the asymptotic information on 
r ni s n , t n coming from Lemma IT^ as well as the asymptotic expansion 

dIT3|) for V (1) ( x), 

one sees that r' ~ 1/(2 logn) and s' ~ 1/(2 logn) as n — * oo. 

It is slightly more delicate to determine the asymptotics of t' , because a cancellation 
of largest terms occurs in the numerator determinant of the formula 



(dE x dE 1 



det 



Or 


ds 


dE 2 


dE 2 


dr 


ds 


dE 3 


0E 3 



•i - <0 (1) (n - 
.l_^,(i)( n _ 

— (n — r 



t + 1) 



t + 1) 



t + 1} 



det J 



(7.17) 



/ 



\ 9r 9s 

for t' that results from Cramer's rule. In this case, in addition to Lemma ITU and (|7.13j) . 
we use the expansion 

1 



(1), 



t+1) 



n — r — s — t + 1 

1 A r s 
- 1+ -+ -+ 

n \ n n 



oil 



n- 



s \ ■ 



- + - 

n n/ 



+ 



1 



x nlog n, 

the first line coming from (|7.13|) . the second resulting from expanding the fraction on 
the right-hand side up to two terms and using the already known fact that t ~ 2 log n. If 
this is substituted in the numerator of (|7.17jl . the claimed estimate t' ~ 2/n is obtained 
after a (tedious) routine calculation. 

For the proof of (|7.15J) . we consider the first and second component in (|7.16jl . The 
first component yields the equation 



1 

r 



nr 



r s 
- + - 

r s 



s 
r 



r's t' 
^r + - 



r't 



+ r> (1) (r + 1) 



r 



n — r — s 



+ r'lfj^ (n — r 



t + l) + sV (1) (n-r 



t + l)+tV (1) (™~ r 



t + 1) 
-t + 1). 



Using the asymptotic information (J7.8|) and (|7.14|) that is already available to us, and, 
in addition, the asymptotic expansion (|7.13jl for the first polygamma function, we infer 
that, as n — > oo, we have 

1 nr' 2 

r r z n 

which, again in view of ([7.8)1 . is equivalent to the first claim in (|7.15jl . The proof of the 
second claim starts with the second component of (|7.16j) . and is otherwise completely 
analogous. □ 
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As an immediate consequence of (|7.14j) and the mean value theorem, we can asymp- 
totically estimate the differences r n — r n _i and s n — s n _i. 

Corollary 17. Let (r n , s n , t n ) G f2 n 6e suc/i that F(r n , s n , t n ) is maximal on fl n . Then 

' ( ' 1 (7.18) 



and 



Sn S n — i 



2 logn 
1 



+ o 



logn 
1 



(7.19) 



2 log n \ log n 

Next, we derive asymptotic estimates for the second derivatives of r n , s n , and t n . 

Lemma 18. Let (r n ,s n ,t n ) G Q n be such that F(r n , s n ,t n ) is maximal on Q n . Then 

1 . ... 2 



r ~ s ~ 



and t'L 



as n — ► oo. 



(7.20) 



2nlog n 

Proof. We proceed in a manner analogous to the proof of Lemma Differentiating 
both sides of ()7.16|) with respect to n, we obtain 




- (1 - r' - s' - t')^ (2) (n-r-s-t+l)\ 
^-(l-r'-s'- t')^ 2 \n -r-s-t + 1) 
\ -(1 -r' - s' - t')4> {2 \n - r - s - t + 1) j 




(7.2i; 



where J' is the derivative of J with respect to n. Using the asymptotic information on 
r n , s n ,t n and r' n , s' n ,t' n from Lemmas ITU and (fTB]) . as well as the asymptotic expansion 

x z \x a 

one finds that the right-hand side of (|7.21|) is asymptotically equal to 



\ n 2 log n f 

Solving the system (|7.21|) for r", s" and t", and using the above estimate, the fact that 
det J ~ 2/ logn from the proof of Lemma HT)| and again the asymptotic information on 
r n ,s n , t n provided by Lemma H3J we obtain the claimed estimates for r", s", and t!' n . □ 

Lemma 19. Given n G N, consider the function <& n (k) of a real variable k ranging 
over 1 < k < n — 1, given by $„(/c) := g{k) + g{n — k), where 

Jv^ \z r v h /c ! 

g(k) = + — + 2r k + 2s k + t k - k\ogs k + - log A; - 21oglogfc. (7.22) 

Sk s k 2s k 2 

Then, for sufficiently large n, <& n (k) has no maxima in the range [l,n — 1], and it 
attains its unique minimum in [l,n — 1] at k = n/2, which is at the same time the 
global minimum in this range. 
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Proof. We begin by gathering information on extremal points of the function $ n (/c). 
Its derivative is <&' n (k) = g'{k) — g'(n — k), where 



^ _ _ T ± _ t[k + kl " kS 'k | 1 ks 'k 

Sk $k ^k 2Sfc ^^k 



+ 2r' k + 2s' k + t' k -logs k -^ + 



Sk 2k k log k 

Using the asymptotic information on r k , Sk, t k and r' k , s' k , t' k provided by Lemmas [T41 and 
HE! respectively, we see that 

g'(k) ~ log A;, as k — > oo. 

Now let fc„ be a solution of the equation &' n (k) = 0. If k n <C n, then gr(n — k n ) 
dominates g(k n ), hence $' n (k n ) = g'{k n ) — g'(n— k n ) cannot be for large n, a con- 
tradiction to the definition of k n . On the other hand, if k n x n, then we infer that 
= g'(k n ) — g'(n — k n ) ~ log k n — log(n — k n ), whence k n ~ n/2. 

Next, we claim that, for large n, any such solution k n is a local minimum. To see this, 
we must compute the second derivative of $> n (k). We have = g"(k) + g"{n — k), 

where 

_ 2 Aks' k k 2 s'l 2k 2 (s' k ) 2 2r' k 2r k s' k kr" k 2kr' k s' k 

9 — ~ 2 ' 3 ' 2 ' 2 

Sk S k S k S k Sk S k Sk s k 

krks" k 2krk(s' k ) 2 s' k ks k k(s' k ) 2 „ n s' k 

i 2 — 3" "2" o — 2~ 3 ft ft ft — 

^ksl k{s' k ) 2 1 2 2 

S fc + S 2 2A; 2 + k 2 \0gk + k 2 (hgk) 2 ' ^' ^ 

Using the asymptotic information on r k , Sk, t k , r' k , s' k , t' k , and r' k , s k , t" k provided by Lem- 
mas O EE and HS1 respectively, we see that 

g"(k) ~ -, as k — > 00. 

In particular, the estimate ()7.15|) has to be used to approximate the terms 

2 Aks' k 2k 2 (s' k ) 2 _ 2 2 

2 + 3 — 3 [Sk kS k ) 

S ft S k S k S k 

occurring in (|7.23j) . Thus, if k n ~ n/2, we have 

= g"(k n ) + g"(n - k n ) ~ 2/A; n ~ 4/n. 

It follows that $"(&;„) > for sufficiently large n, which establishes our claim. 

Now we suppose that there were two sequences (k n ) and (k n ) of local minima of 
$ n (fc), with the property that there is no n Q such that the subsequences (k n ) n > no and 
(k n ) n > no are identical. Then there exists a sequence (ra 3 -)i>o with A; nj 7^ k nj for all j. 
Since both k nj and A; n . are local minima, there must be a local maximum in between, 
m nj say. This gives a sequence (m nj )j> of local maxima of $ n . (fc). Since all these m nj 
must be solutions of & (k) = 0, this contradicts our previous finding that solutions to 
& n (k) = satisfy $"(A:) > for all sufficiently large n. Therefore, there cannot be two 
essentially different sequences of local minima. 
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Trivially, we have $'(n/2) = 0. Hence, for large n, the only local minimum of $(&) is 
k = n/2. Since the above argument also shows that there cannot be a local maximum 
for large n, this local minimum is at the same time the global minimum. □ 

8. Auxiliary lemmas II: Approximation of the sum J2F( r i s >*) 

In the previous section we determined the location of the point (r n , s n , t n ) where the 
function F(r,s,t), which we defined in ()6.3|) . and which is equal to the summand of 
the sum ()4.1|) that we want to approximate, attains its maximum. Furthermore, we 
derived asymptotic properties of r n , s n , and t n . In this section we now embark on 
approximating the sum itself, thereby heavily relying on the information found in the 
previous section. 

We begin by stating (well-known) approximations for sums of exponentials of qua- 
dratic expressions in the summation indices. We provide proofs, which take advantage 
of the Poisson summation formula, for the sake of completeness. 

Lemma 20. (i) Let a, (3 £ R with a > 0. Then we have 

oo 2 / — 

exp(-cm 2 - (3n) = J— (l + 0(a 



as a 



as a, 7 



with an implied constant depending only on a. 

(ii) Let a, /3, 7 £ R with a > and (3 2 — 4<yy < 0. Then we have 
00 00 2 

n=— 00 m=— 00 V ' ' 

with an implied constant depending only on a and 7. 
Proof. We use the Poisson summation formula 

00 00 poo 

Y^ f( n + a )= ^Y, ex P(27r^o) / f(t)exp(—2nifit)dt, (8.1) 

n=— 00 /i=— 00 ~ 00 

which holds in particular for continuous / of bounded variation such that fit) — > as 
t -> 00 and \f(t)\ dt < 00; cf. [221 Thm. 45], (3.12.1)], or [HI (5.75)]. 

(i) Setting f(t) = exp(— at 2 — (3t) and a = in (|8.1j) . and using the Gaufi integral 
evaluation, we obtain 

00 00 ^.00 

exp(— an 2 — /3n) = / exp(— at 2 — f3t — 2irint) dt 

n=— 00 /i=— 00 — 00 

1 £ ex P (< g+ 4 2 ;^ 2 

/-*=— 00 ^ 



where the last line results from isolating the term for fi = from the rest of the sum in 
the next-to-last line. 

(ii) This follows by applying part (i) successively to the summations over n and 
m. □ 
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The next two lemmas are of a technical nature. They are needed in the upcoming 
estimations. 

Lemma 21. For r > \ and T > we have 

9(r ) :=T ( 1 + ^_( T + T+ 0i og ( 1 + ^ <o . 

Moreover, for any fixed t, the function g{T) is monotone decreasing in T . 
Proof. We compute 

T ( T 

g'(T) = — — ^ - log 1 + - 

yV ; 2r(r + T) & V r 
and 

= l-2r-2T 
9 { ' 2(r + T) 2 ' 
Our assumptions on r and T imply g"{T) < for all T > 0. This, together with 
g'(0) = 0, implies g'(T) < for T > 0. Since also g(0) = 0, the latter fact implies 
g{T) < for all T > 0. □ 

Lemma 22. Let r = r n + i? ; s = s n + 5, and t = t n + T . Then we have 

F(r + l,s,t) = ( 1 \ ™^ (n-r-s-t) 
F(r,s,t) \ r) s(r + l) 2 

( R , ^ R+S+T 

_ 21og f 1 + ^±lVlogfl + ^ + logfl- ^ + " + T 



r n / \ s n / \ Ti r n S n t n 

+ ° ( >„ + i<>) 2 + ^ + n-r n -a n -t n -/2-5-T) ) ' (8 ' 2) 
F(r, S + 1, t) ^ + 1 ^ n — s+1 (n-r-s-t) 



F(r,s,t) V S J r(s + l) 2 

f g < ^ + g 

= exp l-^T^) (n - r ^ Sn + 1) -^ 

-i og fi + ^-2io g ri + ^Vio g fi * + * +T 



^ '"n s n 



^ ( s n "I" 'S') 2 s n Tl r n S n t n R S T J ^ 



and 

F(r, s, t + 1) (n — r — s — £) 
F(r,s,t) ~ /•(/ - 1) 



=exp (i" iog ( i+ 9" iog ( i+I ^ 

+ logfl- R+S + T \+o(± + : -]). (8.4) 



ti Tn &n t n J \nn r n s n t n 
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Proof. All formulae result from a straightforward application of Stirling's formula, com- 
bined with Equations (|7.5|) - (|7.7|) and Equation (|7.4|) . □ 

We shall now estimate the contributions to the sum J2 0<rs t<n F(r, s, t) which come 
from triples (r, s,t) that are "far away" from (r„, s n , t n ). The precise statement is as 
follows. 

Lemma 23. Given n e N, let 
Qn '■= Ur, s,t) e M 3 : \r-r n \<\/n, |s - s n \ < \/n, \t - t n \ < \/\og n log log nj. 



F(r,s,t) = - — — F(r n , s„,t n )o(l), as n -> oo. (8.5) 
r s log n 

(r,s,t)en„ 1 " 



Proof. In the subsequent computations we set r = r n + R , s = s n + S, and t = t n + T. 
Using Stirling's formula and Equations (|7.5|) - (|7.7|) . together with the estimate (|7.4|) for 
the digamma function, we find that as n — > oo 

F(r n + R,s n + S, t n + T) = F(r n , s n , t n ) 

( R, 

x exp [n - r n - s n - t n ) 

V T n 

(R\ S 
1 H J (« - r n - s n + 1) 

+ {n - r n - 2s n - t n - R - 2S + §) log ^1 + 



T T 



2,t n 



(t n + T + i)logM 



- (n - r n - s n - t n - R - S - T + ~) log ( 1 'Z~~~~~r I 



1 1 1 



r n + R s n + S t n + T 



i + R + 1L + L 

n r n s n t n R S T r n s n t n 

as long as r, s, t — > oo for n — > oo. 

We shall use this formula mainly for i? <^ r n , S C s„ and T C n, in which case we 
may write 

, , R\ R R 2 /n /i? 3 
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and similarly for log(l + — ) and log(l R±s±r_\ Thus we obtain 

F(r n + R,s n + S, t n + T) = F(r n , s n , t n ) 
IP 
2rl 



x exp ( — -^-^{n — 2r n — s n — t n — 2R — S — T — |) 



n 



S 2 , -i \ R(2R + S + T+£ 

(n - r n - 2s n - t n - R - 2S 1 ^ V 



n 



S{R + 2S-\) {R + S + T)(R + S + T-l 



T 



+r ( i+ i)- (i " +r+ M i+ £)- 



V2tl 



r n + R s n + S t n + T n — r n — s n — t n — R — S — T 
R S T nR 3 nS 3 (R + S + T) 2 

r n s n t* r* si n 



16) 



In order to establish our claim, we have to cut the region Q n \Q n into several pieces, 
which we consider separately. 

First, let R = y/n, \S\ < \Jn and < T < log 2 n. Then it follows from (|8.6J) together 
with Lemmas ITU and I2T1 that 

F{[r n + R\, [s n + S\ 1 [t n + T\) =F (r n + R', s n + S', t n + T) 

(TV^ ( ft 

-g-J + ° ( ~ 

<F{r n ,s n ,t n ) exp(-2(logn) 2 (l + o(l))). (8.7) 
Furthermore, for R = y/n, \S\ < y/n, and T > log 2 n, we have, for large n, that 

< 1, (8.8) 



F(r, s,t + l) n — r — s — t 



F(r,s,t) r(t+l) 

and, thus, by induction on t, it follows that (|8.7|) holds for all T > 0. 

Next we let T < 0, but such that t = t n + T ^> 1. Still assuming R = \Jn and 
\S\ < y/n, we obtain from ()8.6|) again the conclusion (|8.7|) . 

Finally, if i? = \/n, |5| < ^/n, but t <C logn (which now also includes very small t), 
then, by Lemma El 

F(r,s,t+1) n — r — s — t 
F(r,s,t) ~ r(t + 1) > 
for n large enough. Hence, the inequality (|8.7jl will be true for all t = t n +T < log logn, 
say, once it is true for t — t n + T — log log n, which however we already know to be the 
case. 

In summary, up to now we have established (J8.7)) for R = y/n, \S\ < y/n, and all T. 
Now, for R > y/n, \S\ < y/n, and all T, by (|8.2j) and Lemma ITU we have 

F(r + l,s,t) ( 41og 2 n / , NN \ 

; <exp « (i + (i)) < 1. 



F(r,s,t) V v 7 ™ 
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Thus, by induction on r, we see that (|8.7|) holds for R > y/n, \S\ < ^/n, and all T. 
Since there at most 0(n 3 ) such triples (R, S,T), we obtain 

£ F(r,s,t) = o(F(r n ,s n ,t n )). (8.9) 

r-r n >\fn 
\s—s n \<y/ri 
0<t<n 

A similar inductive argument using ()8.2|) and ()8.3|) shows that 

F(r n + R,s n + S, t n + T) < F(r n + y/n, s n + yfn, t) 

for R > s/n and S > y/n, and hence, because (j8.7|) holds for R = S = y/n and arbitrary 
T, that ()8.7j) also holds for R > y/n, S > y/n, and all T. Again, there are at most 
C(n 3 ) such triples (R,S,T), and thus 

£ F(r,s,t) =o(F(r n ,s n ,t n )). (8.10) 

0<t<n 

A completely analogous argument shows that ()8.7|) also holds for R = —\fn, \S\ < 
■s/n, and arbitrary T. Now, for R < —\/n, \S\ < y/n, and T < ^ l /n 1 by ()8.2|) and 
Lemma ITU we have 

F{r + l,s,t) f 4\og 2 n n \ 

> exp -=— (1 + o(l)J > 1. 



F(r,s,t) J V v 7 ™ 

Thus, by a reverse induction on r, we see that ()8.7|) holds for i? < —y/n, \S\ < \fri and 
T < yfn. If, on the other hand, R < —\fn, \S\ < y/n, but T > y/n, then in ()8.6|) the 
term — T/12t^ is of asymptotic order > \fnj log 2 n. By repeating previous arguments, 
it then follows that 

F(r n + R,s n + S, t n + T) < F(r n , s n , t n ) exp f -— ^1 + o(l)) 



. 12^ 

< F(r n , s n , t n ) exp (-2 log 2 n(l + o(l))) . 

That is, inequality ()8.7j) holds for R < —\/n, \S\ < \fn, and all T. Since there are at 
most C(n 3 ) such triples (R, S,T), we obtain 

Hr, s, t) = o{F{r n , s n , t n )). (8.11) 

r—r n <—*Jn 
\s—s„\<y/n 
0<t<n 

A similar argument using ()8.2|) and (j8.3|) shows that 

F(r n + R,s n + S, t n + T) < F(r n - y/n, s n - y/n, t) 

for R < —y/n and S < —y/n, and hence, because (j8.7J) holds for R = S = —\fn and 
arbitrary T, that ()8.7|) also holds for R < — -y/n, S* < —y/n, and all T. Thus, 

£ F(r,s,t) = o(F(r nj s n ,t n )). (8.12) 

r—r n <—*fn 
s—s n <—y/n 
0<t<n 
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The term F(r,s,t) is almost symmetric in r and s. It is not difficult to convince 
oneself that minor modifications in the above arguments lead to the conclusion that 

F(r,s,t)+ F ( r '^)+ Yl F(r,s,t) = o(F(r n ,s n ,t n )). (8.13) 

\r-r n \<y/n r-r n <-^/n r-r n >^/n 

\s—s n \>^/n s—s n >^/n s—s n < — ^/n 

0<t<n 0<t<n 0<t<n 

The range which remains to be considered is \R\ < y/n, \S\ < s/ri, \T\ > ^/{ogn\og\ogn. 
From ()8.6|) we obtain 



F(r n + R,s n + S, t n ± ^Jlogn log log n) 

R 2 
2r 



< F(r n , s n , t n ) exp ( - -^{n - 2r n - s n - t % 



^( n - r R(2R+S) S(R + 2S) logn(loglogn) 2 



n 



(8.14) 



Furthermore, if \R\ < y/n, \S\ < y^n, and T > ^lognlog logn, we have 



F(r,s,t + 1) ( Vlognloglogn , 
F(r,s t t) < 6XP { T n (1 + » ' (8 ' 15) 



and if T < — i/lognloglogn, 



F(r,s,t+1) ( Vlognlog logra . 

> exp (1 + o(l)) . 



F(r,s,t) V U 

Hence, writing q for the right-hand side of (|8.15jh we see that 

F(r,s,t) <F(r,s,t n + v/loi^ log log g^-v^iogiogn 



for t > t n + y/\ogn log log n and 

F(r,s,t) < F (r,s,t n - ^log n log log nj q ^-V^io g \ 



; n—t 



for t < t n — \/logn log logra. In combination with (|8.14j) this implies that, for fixed r 
and s with \R\ < y/n and |5| < \/n, we have 

\t-t n \>y/\ogn\og logn 

/ ^>2 

< F(r n , s n ,t n ) expf - — (n - 2r n - s n - t n ) 

\ n 

S 2 , s R(2R + S) S(R + 2S) „. A 
- — (n - r n - 2s n ) - 1 - -± + 0(1) 

loenftoeloEtnl 
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Using the definition of q to approximate the term 2/(1 — q), we find that the last line 
in the above expression is 

/ (log log nf \ 
6XP { 2 ^ °^ ^ J ' 

Therefore, we obtain 

|r— rnl^y'n 
|t-tn|>v / I°grilog logn 



<F(r n ,s n ,t n ) ^ exp( 

0<r,s<n V 



p2 c2 



i?(2i? + 5) + IS) 



Applying Lemma EH (ii), and using again Lemma El we obtain 



oil). 



V F(r, s, t) < F{r n , s n , t n )-^o(l). (8.16) 

\s-s n \<^n 



\t—tn\> Vlog™ l°g l°g n 

Finally, if we add up (fOjl . (I8~TU1) . (f8~TTTl. (l8~T2l . (jHH, and (fFTTKj) . the claim follows. 

□ 

The benefit of Lemma ESI is that 

J2 F{r,s,t)= Y] F{r,s,t) + F{r n ,s n ,t n )-^-o{l). (8.17) 
* — z — log n 

0<r,s,t<n (r,s,t)e<3„ 6 

As we shall see presently, the first term on the right-hand side dominates the second. 
More precisely, we now prove the following statement. 

Theorem 24. We have 

l \ t-i/ i \ (2tt) ^ r n s n\/tn 

c(n) ~ F(r n , s„,t n ) , as n -> oo, 

n 

where (r n ,s n ,t n ) is the unique solution of the system of equations (|7.5j) - (|7.7j) . 

Proof. If (r, s,t) G Q n , then we deduce from f!8.6|) that 

/ R 2 ( fr 2 
F(r, s, t) = F(r n , s n , t n ) exp f - — + 2r n - s n - t n + ^ 



2 2 \ IL • It' I II I 1 1 — 1 

|!(i+ (i))+ (i)Y 
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As long as (r,s,t) G Q n , we have for |e r |, |e g |, \e t \ < 1 that 
F(r+e r , s + e s ,t + e t ) 

= F(r n , s n , t n ) exp ^ - — + 2r n - s n - t n + O 



2si \ \ n 



(R + e r )(S + £s )[- + - + 0[-)) - (r j" £t)2 (l + o(l)) + o(l) 



nil 2t r 



n J n 

R 2 ( fr 2 
— \n + 2r n -s n -t n + 0\- 



F{r n , s n , t n ) exp ( - — ( n + 2r n - s n - t n + O [ — 



^(n-r n + 2s n + o(^])-Rs(- + - + o(- 



_rj|-- it, • ill i — I II — 1 ■ ■ — I 

2s z n \ \n J J \r n s n \ n 

T 2 



^■(l + o(l))J (l + o(l)). 

Thus, on the right-hand side of (J8.17|) . instead of summing over (r, s, t) G Q n , we may 
sum over (R,S,T) G Z 3 with \R\ < y/n, \S\ < y/n, \T\ < ylognloglogn. Therefore, 
the first term on the right-hand side of (|8.17|) is equal to 

/ R 2 / f r 2 
F(r n ,s n ,t n ) ^2 expi - — ln + 2r n - s n -t n + ( — 

\R\<y/n V r n V V n 

\s\<vn 

|T|<vTognlog logra 

S ~ f n-r n + 2s n + o( S ^-X\-Rs(- + — + O n 



2s l n \ \ n J J \r n s n \n 

+ + (8-18) 

Next we argue that we may extend this sum to all (R, S, T) G Z 3 without introducing 
significant errors. More precisely, we claim that 



^ eXp ( " || ( n + 2r " ~ Sn ~ £ " + ( n )) 



V f n - r n + 2s n + O ( ^X] - RS (— + — + O ( ' 



2n I ' ~ ■ il • / £■ ■ — i if I 1 — i 

Si V \ n J J \ r n S n \Tl 



^(1 + 0(1))) (1 + o(l)) = o (-\- ) , (8.19) 
2t n J \log n 



where Yl is over au •S^) e ^ 3 which violate at least one of \R\ < y/n, \S\ < \fn 
and \T\ < \/log n log log n. This can be proved by using arguments similar to but 
significantly simpler than those in the proof of Lemma 1221 We leave the details to the 
reader. 
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Finally, we collect our findings. By the definition of c(n), (|8.17|) . (|8.18|) . and (|8.19|) . 

we have 

c ( n ) = F ( r ' s '*) 

0<r,s,t<n 

= F(r n ,s n ,t n ) Y ex pf- (n + 2r n - s n -t n + O (—X) 



-oa<R,S,T<oo 

S " f n-r n + 2s n + o( S ^X\ - Rs(- + - + O fl 



1 77 

+ F(r n ,s n ,t n )- — n - o(l) 



Sn V \n J J \T n S n \ Tl 

^(1 + 0(1))) (1 + 0(1)) 



log 2 n 

Now we apply Lemma EZDT i) to the sum over T and Lemma OUIf ii^ to the sums over i? 
and S. Little simplification then yields that 

c(n) = F(r B , s„, g^^ (1 + (i)) + F ^ Snj tn )™ {l). 

n log n 

By the asymptotic information on r n , s n , t n provided by Lemma IT^ we see that the first 
term on the right-hand side dominates the second, whence the claimed result. □ 

By the preceding theorem, we have now expressed the asymptotics of c(n) in terms 
of the solutions r n ,s n ,t n of the system of equations (|7.5J1 - (|7.7|) and the value of the 
function F(r, s, t) at (r, s, t) = (r n , s n , t n ). While the asymptotics of r n , s n , t n is already 
known from Lemma fill the missing step in the proof of Theorem |H] is to find an explicit 
expression of the asymptotics of F(r n , s n ,t n ). This is done in the next proposition. 

Proposition 25. We have 

F(r n , s n , t n ) = exp ^n(log n - log s n ) + ^- - r n + - 3n + 2r n + 2s n + t n 

1 (2nfr n t n \ 

log h 1 + oil) as n — > oo. 

2 s n / 

Proof. We apply Stirling's formula to the defining expression for F(r n , s n ,t n ) to obtain 
F(r n , s n , tn) = exp( ( n + \ J log n + (n - r n - s n - t n ) log r n 



■■■in - r n - s n + 1) log s n - yr n + -j log r n - [ s„ + - ) log >>, 
l\. / 1' 



t „ + ^logt„-^- r „- S „- t „ + -|.o g(n - r „- S „- ( „) 
ilog(2jf + o(l)Y 
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Now we make use of the Equations (|7.5|) - (|7.7|) . of Equation (|7.4j) for the digamma 
function, and of the expansion log(l + x) = x + 0(x 2 ). This gives 

F(r n ,s n ,t n ) = exp ( - n + r n + s n + t n + O ( — ) + O ( Tn + Sn — — ) 

\2r n \r n J \n- r n - s n -t n J 



+ -^-n + r n + s n -l + - j+T^T + O r 



t II 



n(n — r n — s n + 1) , n 1 (r n 

+ - " ! --n\ogs n - — + - + 

1 , (2nfr n t n . . 

+ n\ogn- - log h o(l) 

z s n 

The result follows now upon appealing to the asymptotic information on r n , s n ,t n pro- 
vided by Lemma ITU □ 
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